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SMALL  DEFORMATIONS  OF  A PLASTIC- 
RIGID  BODY  AT  THE  YIELD  POINT 
By  E.  T.  Onat* 

1,  Introduction,  Consider  a body  of  plastic-rigid  material 
which  is  loaded  by  nonotonically  increased  surface  tractions. 

Let  us  assume  that  the  history  of  loading  is  completely  specified. 
At  the  earlier  stages  of  loading  the  body  will  generally  remain 
rigid.  If  the  loads  are  increased  further  some  portions  of  the 
body  may  become  locally  plastic.  However,  plastic  deformation 
of  these  portions  may  be  prevented  by  adjacent  rigid  portions 
until  the  moment  when  the  body  actually  begins  to  deform.  Follow- 
ing Hill  [1]^,  we  refer  to  this  critical  moment  as  the  yield 
point  of  the  plastic-rigid  body.  If  the  material  is  perfectly 
plastic  (i.e. , non- strain- hardening)  and  if  the  accompanying 
change  in  geometry  is  disregarded,  plastic  flew  is  found  to  con- 
tinue under  constant  loads  [2]«  However,  when  the  change  in 
geometry  and  the  effects  of  strain-hardening  are  taken  into  ac- 
count, continuing  deformation  under  constant  external  forces  is 
possible  only  in  exceptional  cases  [l  3.  Thus,  as  a rule,  quasi- 
static flow  requires  either  increasing  or  decreasing  external 
forces.  Under  the  assumed  monotonically  increasing  loading,  the 
first  case  represents  a stable  process,  whereas  the  second  case 
leads  to  sudden  collapse. 

* F.  3.  Jewett  Fellow  at  Brown  University,  Providence,  R.  I. 

# Numbers  in  square  brackets  refer  to  the  Bibliography  at  the 
end  of  the  paper. 
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The  distribution  of  stress  and  the  corresponding  velocity 
field  at  the  yield-point  can  be  found  in  some  cases  (e.g, , by  the 
method  of  characteristics  in  plane  strain  problems)  and  in  other 
cases  approximations  from  above  and  below  to  the  yield-ncint  load 
can  be  obtained  [1,2].  However,  in  many  cases  the  present  first 
order  theory  does  not  furnish  a clear  criterion  for  collapse. 
Therefore  it  seems  worthwhile  to  investigate  the  second  order 
effects  in  the  deformation  of  a plastic-rigid  body  at  the  yield 
point. 

In  the  following  the  investigation  is  limited  to  problems 
where  the  entire  body  is  in  the  plastic  state  at  the  yield  point 
and  where  the  incipient  stress  and  velocity  fields  are  known.  A 
perturbation  scheme  is  developed  to  obtain  equations  governing 
Lagrangian  stress  and  velocity  increments  of  the  subsequent  quasi- 
static motion.  Some  examples  are  worked  out  (compression  of  block 
between  rough  plates  and  the  thick  walled  tube  under  internal 
pressure).  In  the  first  case  the  results  obtained  with  the 
theory  are  found  to  agree  with  the  experimental  results.  In  the 
second  example  the  minimum  rate  of  strain-hardening  necessary  to 
prevent  collapse  is  found, 

2*  Notations  and  Stress- strain .Relations  Employed  in  the  Analysis. 
Let  xj  denote  the  rectangular  Cartesian  coordinates  of  a generic 
particle  at  the  time  t.  If  the  velocity  of  this  particle  is 
denoted  by  vj_,  the  velocity  strain  tensor  is  given  by 


where  Latin  indices  take  the  values  1,  2,  3. 
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The  stress  and  deviator lc  stress  tensors  are  denoted 
respectively  by  d^  and  s^: 

siJ  = 8ij  * (3  <’kk)(,i3>  (2) 

where  summation  convention  is  used  and  6^  represents  the 
Kronecker  delta# 

If  the  second  invariant  of  the  deviatoric  stress  tensor 
is  denoted  by  J2, 

J2  = 2 SlJSiJ*  (3) 

The  particular  stress-strain  relations  considered  in  this  paper 


are  DJ_ 

eij  = P(J2)siJ  Dt  whea 


2 DJ2 

J2  2 k and  __  £ 0 


when 


J2  < k2  or  ££2  < C 


Dt 

DJ, 

Dt 


L (*f) 


where  DJ2/Dt  is  the  material  derivative  of  J2* 
The  limiting  case  where 

DJd 

P(j23_»co  ^-^0 

defines  an  ideally  plastic  material; 

Jo  = k2, 


13 

= Xsu 

when 

\ 2 0 

13 

= 0 

when 

\ < 0, 

1 

!- 


(5) 
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The  stress  and  deviatoric  stress  tensors  are  denoted 
respectively  by  and  s^: 


1J  = °ij  ' (3  <kk)  61J  ’ 


(2) 


where  summation  convention  is  used  and  5^  represents  the 
Kronecker  delta# 

If  the  second  invariant  of  the  deviatoric  stress  tensor 
is  denoted  by  J2, 

J2  ~ 2 SijSij* 


The  particular  stress-strain  relations  considered  in  this  paper 


are 

DJo 

eij  = F(J2)sij  Dt^ 

when 

J2 

dj2 

and  — 

dF 

c 

O 

A| 

. (W 

eij  = ° 

when 

J2  < k2 

dj9 

or  —Jk 
Dt 

1 

< 0 

o 

where  DJ2/Dt  is  the  material  derivative  of  J2* 
The  limiting  case  where 

DJ0 

F(J2>— >oo 

defines  an  ideally  plastic  material: 


= Xsij 

when 

X 

2°  ] 

> (5) 

U 

= 0 

when 

X 

< 0. 

Vj 
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3*  Formulation  of  the  Problem.  Consider  a plastic-rigid  body 
at  the  yield  point.  Suppose  that  the  boundary  conditions  are 
such  that  every  element  of  the  material  is  in  the  plastic  state* 
n’oreover  assume  that  a velocity  field  v°  is  associated  with  the 
stress-state  where  and  are  the  Cartesian  components 
of  the  velocity  and  stress  tensors*  v°  and  satisfy  the  bound 
ary  conditions  imposed  on  the  displacements  and  the  stresses  in 
addition  to  the  following  equations: 


S ^ «5  ® — O If  ^ 

sijsir 


v°  ) = 


e = \s° 
ij  U 


where 


X > 0. 


(Equilibrium) 

(Yield  condition) 

(Stress  strain  rela- 
tion) 


(6) 

(7) 

(8) 


We  note  that  this  solution  is  valid  for  an  ideally  plas- 
tic material  and  also  for  a newly  annealed  rigid-strain-hardening 
material  which  has  the  same  yield  stress,  k,  in  pure  shear* 

Wow  assume  that  the  surface  tractions  are  varied  in  such 
a manner  that  the  further  distortion  following  the  reaching  of 
the  yield  point  takes  place  in  a quasi-static  manner. 

We  now  consider  the  early  stages  of  the  distortion  such 
that  the  state  of  velocity-strain  and  stress  differs  only  slightly 
from  the  state  which  is  given  by  Eqs*  (6),  (7),  and  (8). 

Let  be  the  Cartesian  coordinates  of  the  material  point 


ap  at  time  t: 


^ = xA[a  , 6(t)  ] 


(9) 


where  6(t)  is  a monotonously  increasing  function  of  time  and 
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6(0)  =0,  xja,  , 0 1 = a<  • 

If  the  velocity  and  acceleration  at  time  rero  are  con- 
tinuous functions  of  the  space  variables,  the  validity  of  the 
following  expansion  may  be  assumed: 

2 

*i  = ai  + v°(ap)&  + Yi^ap^  V + •••  (10) 

where  to  obtain  a quasi-static  flow,  d6/dt;  d 6/dt  are  supnosed 
to  be  as  small  as  desired.  Since  the  stress-strain  relations 
envisaged  in  this  theory  are  homogeneous  in  time,  the  actual  form 
of  these  functions  is  of  no  importance  in  the  problem  considered. 
The  velocity  of  the  material  point  ap  at  time  6 referred 
to  the  fixed  rectangular  Cartesian  coordinates  is  then  given  by 

= (vj  + + •••)  jr  . (11) 


is 


The  velocity  strain  tensor  defined  in  the  Eulerian  manner 


- + Hi  1% 

' 2^ 9a  , i + 8a_  •" 

p ax  Pi 


(12) 


From  (10),  9ap/9xj  can  be  evaluated  (see  Appendix)  in 
terms  of  v^: 

U3) 


9a  n o 

= 6 r , - 6v°  . + . . . 
9x^  PJ  P,J 


where  a comma  denotes  partial  differentiation  with  respect  to  a^, 
From  the  last  three  equations,  we  obtain 


r , I 1 d6 

eij  = [eiJ  + 6eij  + ] d^ 


(14) 


where  = ^(vi,J  + vj,i^  and 
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'll  = i(U,3  + YJ,l)  - ^(vl|PvP)3  + v3)PvPil)‘  (15> 

Let  {j^j(ap,6)  be  the  components  of  the  stress  tensor  at 

the  point  x^  = x^(a  ,6),  referred  to  the  fixed  rectangular  Car- 

P 

tesian  coordinates.  For  small  deviations  from  the  yield  point 
state  we  again  assume  that 


d,  < = d?.,(a  ) + bo’  .(a  ) + ...  • (16) 

ij  ij  P ij  p 

We  now  examine  conditions  and  equations  which  must  be  satisfied 
during  the  quasi-static  distortion  of  the  body  after  the  yield 
point* 


a.  Equations  of  Equilibrium;  Since  a quasi- static  notion 
of  the  body  is  considered,  the  Equations  of  Equilibrium  must  be 
satisfied  during  the  subsequent  flow: 


9o . . 0d,  . da„ 

— ±i  = — U . — S s 0, 

ax3  eaP  ax3 

where  body  forces  are  neglected. 

Using  (13)  and  (16)  this  becomes 


(17) 


ili  ♦ 6[iy 

8ej  aSj 


- V 


n 


a°?3 


P,J  9a 


] = 0 


(18) 


where  only  the  terms  up  to  the  order  b are  retained.  The  first 
term  of  this  equation  vanishes  according  to  (6),  hence  the  Equili- 
brium Equations  satisfied  by  the  Lagrangian  stress  increments  are 

^3  . v°  aa?3  . 0 


(19) 
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b.  Stress-Strain  Relations.  An  investigation  of  the  stress 
strain  relations  at  time  zero  and  time  b provides  six  additional 
equations  for  Lagrangian  Stress-  tixid  Strain  increments* 

At  time  6=0, 


where 


dJ5(a  ,6) 
■°  = F(J2)s°  - 2 

J-J  d ij  d b 


+ 26s. .3, 


2'“p,w'  - ij-ij 


(20) 


and  denote  increments  of  the  deviatoric  stress  tensor* 
fore  if  the  terms  up  to  the  order  of  unity  are  retained, 


There 


and  (20)  gives 


dJ2  o * 

vr  = Vpq 


O I 

s s 

pq  pq 


*cr 


3iJF(J2) 


(21) 


where  the  right  hand  side  is  not  summed  and  any  combination  of 
(i,j)  may  be  chosen  provided  that  s^  j*  0* 

If  the  material  is  ideally  plastic  (21)  becomes 


s 


pq  pq 


= o. 


(22) 


At  time  6,  in  addition  to  the  equation  of  incompressibility,  the 
stress-strain  relations  provide  four  independent  equations  of 
the  following  type: 


elj  _ £pq 


(23) 


pq 

where  again  the  summation  convention  does  not  apply.  Using  (8), 
(15)  and  the  definition  of  a^  , and  retaining  terms  up  to  the 
order  of  b one  obtains  the  following  equations 
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* 0 . o » o » 

t..S  + 6,  , S „ “ e^n  s.  . 

ij  pq  ij  M P3 


- els?,  = 0, 


:pq  ij 


(24) 


where  (i,j,p,q)  are  to  be  chosen  in  order  to  give  four  independ- 
ent equations. 

On  the  other  hand  the  equation  of  incompressibility,  which 
is  obtained  from  the  contraction  of  the  tensorial  equation  (4), 
gives 


and  using  (14)  and  (15),  and  observing  that  e°^  = 0,  we  get 

(25) 


o o _ 
Y - V.  V , = 0. 
Ti,i  i,P  p,i 


We  note  that  Eqs.  (19),  (21)  (or  22),  (24)  and  (25)  provide  nine 
linear  equations  for  the  nine  unknowns  and  Y^»  Ii  the  case 
of  plane  strain,  Eqs.  (19)  and  (21)  or  (22)  provide  three  inde- 
pendent equations  for  three  unknown  stress-increments.  Therefore 
if  the  boundary  conditions  are  of  the  stress  type  the  problem  is 
statically  determined. 


c.  Boundary  Conditions.  We  first  examine  conditions  which 
must  be  satisfied  at  a boundary  surface,  Sp,  of  the  body  at  which 
surface  forces  are  prescribed.  Denote  by  the  components  of 
the  surface  tractions  in  the  fixed  rectangular  Cartesian  coordi- 
nate system: 

QjUpt&l  = Q|<a  ) + 6Qj(a  ) + ...  . (26) 

Suppose  that  Sp  is  defined  by  the  following  equation 

ai»a2  »a3 ) = 0 (27) 


where  as  before 
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or 


xi  = xi<ap,6) 
*p  = ap(x1(6) 


and  it  is  assumed  that  points  originally  on  the  boundary  surface 
remain  on  a surface  during  the  deformation.  At  time  6^,  the 
unit  normal  np  of  the  surface  Sp  has  direction  cosines  propor- 
tional to  9P/6xp.  Using  the  series  expansions  and  (15)  one  ob- 
tains: 

n = n°  + 6n'  + ...  = k 2JL  = k ££-(&  - 6v°  ), 

P P P 9x  8«l  rp  r,p 


or 


% 


0ay 
r 

+ 6n'  = k(M 6F  v°  ). 


9a 


♦ * r,p' 


(28) 


The  magnitude  of  k must  be  such  that 


n n = 1 
P P 


(unit  normal). 


(29) 


From  (13)  and  (14) 


k = + 


F F v°  n 

s 1 + 6 .r>P 


(30) 


where  the  sign  will  be  properly  chosen  to  obtain  the  outward 
normal  of  the  surface. 

On  the  other  hand,  the  surface  tractions  are  related  to 
the  stresses  bv 

^ = VPi 

or  1 f t 

+ + •••  = (dp^  + k^pi  + ...)(n^  + 6n  + ...), 


1 pi  p pi  p 


From  (28)  and  (30); 
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jJP 


)/f  7f 


• J »J 


^ [ ' F>,k 


F P v 

■>i,.».r..M  - F „v°  1 

>r  ‘->P 


and  therefore  the  Increments  of  the  surface  tractions  are  related 
to  the  increments  of  the  stresses  by  the  following  formula: 


Q'  = [f  [LaLi^Lsl  . r 

1 Pi  i/P  .F  . \/iT  F ,P  F ,F  , ’ r,P 

v ,j  ,j  i J »j 


VP,3F,1 


- F.?V°  „].  (31) 


As  to  the  displacement  boundary  conditions,  we  infer  that 
on  Sv,  Yi  is  given  as  functions  of  the  original  coordinates* 

The  theory  has  now  been  taken  to  a point  from  which 
applications  may  be  made  to  specific  problems. 


Example  1.  Compression  of  Block  Between  Rough  Plates: 

Consider  a block  of  strain-hardening  (or  ideally  plastic )- 
rigid  material  compressed  between  rough  plates  (Fig.  1). 

When  the  block  is  very  wide  compared  with  its  height,  the 
slip  line  field  corresponding  to  the  incipient  motion  is  independ- 
ent of  x,  and  the  incipient  stress  and  velocity  fields  are  given 
by  Prandtl  [3]  and  Nadai  [4]: 

n 


co 

k 

£ 


T° 

xy 


= - e - 


+ 2V1 

h 


mx 

h 


- m 


my 

h 


£ 

h2 


(32) 


and 


All- 90 


r\ 


11 


£ 

U 


= c 


v 

V 


* 


m 


2 

L- 

h2 


(33) 


where  k is  the  yield  stress  in  pure  shear,  m is  a measure  of 
roughness  [0  < m < 1;  m = 1 for  perfectly  rough  plates]  and 


-1 


c = -SlO — SB  +\ji  - 


m 


m 


1 


= . L + jalaZia  + 1 \fi- 

h 2 mV 

m 


2 

m • 


(3*) 


These  equations  state  that  when  the  normal  pressure,  p, 
applied  by  the  plates  reaches  the  value  given  by 

£ = c + 32S  + 2 \J  1 - m2 

k h v 

the  plastic  rigid  block  begins  to  deform  according  to  the  velocity 
field  ( 33 ) • As  stated  before  we  are  interested  in  the  subsequent 
neighboring  deformation  of  the  body  and  in  the  rate  of  increase 
of  the  external  forces  to  achieve  this  deformation* 

If  (x,y)  are  regarded  as  Lagrangian  coordinates  and  if 
the  current  coordinates  are  denoted  by  (x,y), 

x = x + 6^  + ••• 

y - y + &■§  + ... 

where  6 denotes  the  compression. 

Equations  (19)  and  (21)  provide  three  independent  rela- 
tions to  determine  the  Lagrangian  stress  increments  dy,  T^y. 
#rora  (19),  (32)  and  (33) 
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5d  1 

dx  1 

2Jm 

\ 

2 

+ — Z2  + 

■°  i 

l 

dx 

ay 

h2 

v (35) 

3d  ' 

dx  * 

y 

3y 

yx 

+ dx 

17  J 

where 

/ — 

— 

2 

f (y) 

= 2 V1  • 

- m2 

y 

From  (21) 

i i 

d - d„ 

x y 

+ iff 

= - 

f.y 

m2k2F 

h . (36) 

■ y 

If  it  is  assumed  that  (dx  - cJy)  and  x^.  are  functions  of  y only, 

the  above  equations,  with  the  boundary  conditions  x = 0 when 

xy 

y = 0,  and  x^  = mk  when  y = h,  give 


x =0 
xy 


dy=|f-2k-^x+b  ^ 


°x  =d--  - 


h 

f*y  h 


i 


(37) 


y m2k2F  y 

where  b is  an  undetermined  constant.  To  determine  b let  us  con- 
sider the  surface  AB  which  after  a small  deformation  becomes  a 
plane,  A'B',  parallel  to  the  y-axis,  Fig.  1. 

Transformation  from  AB  to  A!3’  is  given  by 


£ + r (C  + £-£)  = 2 = constant 
h h h m h 


_ y 

" h * 


(38) 


V 


Since  the  left  hand  edge  of  the  block  is  free  from  stress,  the 
overall  equilibrium  condition  demands  that 

flhd-b/h) 

oxdy  = - mkx.  (39) 

uO 
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Using  (32),  (37),  (38)  the  last  equation  gives 

h 


mkC  + bh  - 


»2k2F 


LlL  dy  = 0, 

y 


hence 


mkC  + bh  = - 


JO 

-2—  arc  sin  ra. 


mk^F 


(40) 


The  normal  stress  in  y-direction  is  then  obtained  from  (37)  and 


+ bo  1 = (f-c)k  - nk  { [1  + r]  - 
v h h h 


(40): 

* = rt°  + Sri ' = (f-c)k  - mk  , 

y h 

The  average  normal  pressure  required  to  continue  the  deformation 

is  therefore  

p = - (2  \/l  - m2  - c)k  + [1  + i]  + — yS-Slfl-B.  (42) 

2h  h h mk2F 

In  order  to  compare  this  result  with  experiment  [5]»  let 


us  consider  the  case  m = 1: 


£ = Z + -k  + k [ _2L_  + -k.  -I 
k 2 2H  h L k3p  2H  J # 


(43) 


The  dimensionless  quantity  Fk^  can  be  obtained  from  the  stress- 
strain  curve  in  compression: 


Fk" 


<£> 


where  o is  the  yield  stress  in  compression  and  do/de  is  the  slope 
of  the  stress-strain  curve. 

The  stress-strain  curve  of  the  material  used  in  the  ex- 
periments is  shown  in  the  Fig.  2.  The  material  begins  to  deform 
plastically  at  50  Bars  and  at  this  point 

-2-  = Its-  ag/df.^ii2?. 

Fk3  1/3  0 
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Hence  according  to  (43) » [the  slope  of  the  load-deforma- 
tion curve ] at  the  beginning  of  the  deformation  is  strongly 
dependent  on  the  rate  of  strain-hardening  since 


The  results  obtained  from  the  experiments  agree  with  this 
observation  (Fig,  3),  Moreover  the  straight  line  (I)  defined  by 
the  equation  (43)  follows  closely  the  experimental  results* 

In  order  to  have  an  idea  on  the  shape  of  load  deformation 
curve,  one  can  approximate  the  stress-strain  relation  with  an 
ideally  plastic  material  (F— ^00)  as  shown  in  Fig*  2.  The 
straight  lines  (II)  defined  by  (43)  are  approximate  load-deforma- 
tion curves  corresponding  to  this  idealized  material  (Fig*  3). 


Mamnis-S*  .OzUnlej:. uq^TL-iat ■Pr.qgsmrs • 

A tube  made  of  plastic-rigid  material  is  subjected  to 
internal  pressure.  The  tube  is  supposed  to  be  long  enough  so 
that  the  axial  deformations  of  the  tube  can  be  neglected*  If 
the  yield  stress  of  the  material  is  denoted  by  k,  incipient  stress 
and  velocity  fields  are  given  by  (Fig,  4) 


and 


= 2k  log  X 
Og  = 2k  (1  + log  X) 


n 


> 


o 

✓(r) 


V 


(44) 


The  material  points  situated  on  the  circle  of  the  radius  r would 
move  after  time  6 to  another  circle  with  the  radius  r 
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r = r + •&&  + .«.  • 
r 

The  physical  comDonents  of  the  stress  in  these  elements 
would  be  given  by 

d = o®  + 6d  +••• 
r r r 

= c©  + + • • • 


Tr© 


= 0. 


The  Equilibrium  Equation  satisfied  by  the  Lagrangian  Stress- 
increments  is  found  either  from  the  tensorial  Equation  (21); 
or  from  a straight  forward  derivation: 

ao*  o'  “ o' 

_£  + -I 2 + 4k  -%  = 0.  (45) 

9r  r r3 

On  the  other  hand,  Equation  (16)  gives 


Eliminating  from  these  Equations 

Oh. & + hte  = 0 

Or  k2Fr3  r3 

■“  + ^ [ 4 - -i-  ] = 0 

ar  r3  k3F 

o'  - ( 4 i-  ] + c = 0. 

2r^  k3F 

Using  the  boundary  condition 


(46) 


at 


r = b 


* This  equation  must  first  be  rewritten  for  curvilinear  coor- 
dinates# 
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we  get 


and 


8r  =f  » -+  1[J5 

r * k3F  " v 


<5  (r,6)  = 2k  log  £ + [4 i-]  [ -L 

r b 2 k3p  r2 

It  is  easily  seen  from  this  equation  that  if 


(47) 


F > 

4k3 


(48) 


I or(a,6)  | < | <5^ ( a , 0 ) | , 

hence  a quasi-static  motion  can  only  be  maintained  with  decreas- 
ing pressure.  [In  the  case  where  F — >oo  (ideally  plastic  mater- 
ial), (47)  checks  with  a previous  analysis  [6]].  However  if 

F < — 

4k3 


deformation  can  be  continued  with  increasing  pressures.  In 
terms  of  a simple  stress-strain  relation  in  compression  the  in- 
equality (48)  can  be  written 

At.  »/o) 

< V3  collapse 

QE 


d(c/oQ) 

de 


> V/3 


no  collapse. 


In  Fig.  (3),  the  limiting  curve 
shown. 


separating  the  two  regions  is 
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Appendix:  It  is  well  known  from  Tensor  Analysis  [?1  that 


ia!  = gpm  asi  (49) 

0x  8am 

where  gpm  are  the  contravariant  components  of  the  metric  tensor 
of  the  curvilinear  system  defined  by  ap  in  the  deformed  position* 
However  g^m  are  given  [8]  to  the  desired  order  of  approximation 

^ gP“  = 6 - 6(v°  + v°  ) 

pm  P»m  m,p 


and  since 


( 4-9 ) gives 


-22L-  = 6 + v?  6 + , 

_ m im  i.m  1 

oa 


= V " iv 


9x 


4 t • • 
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